A  Computer  Model 
of  the  Gypsy  Moth 
and  its  Fungal 
Pathogen 


cokA 

s 

H3 


BY  RONALD  M.  WESELOH 


Abstract 

A  one-host  generation  computer  model,  with  a  time-step  of  1  day,  has  been  developed 
to  simulate  the  interaction  between  the  very  effective  fungal  pathogen,  Entomophaga 
maimaiga,  and  its  host,  the  gypsy  moth.  This  publication  describes  the  structure  of  the 
validated  model,  the  kinds  of  model  inputs  needed  and  the  outputs  produced,  and  gives 
examples  of  model  results.  The  model  can  be  used  to  construct  scenarios  of  future  gypsy 
moth-fungus  interactions  and  to  evaluate  the  effectiveness  of  the  fimgus  as  a  biological 
control  agent. 
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The  gypsy  moth  fungal  pathogen,  Entomophaga 
maimaiga,  was  first  found  to  be  established  in  North 
America  in  1989  (Andreadis  and  Weseloh  1990),  and 
positively  identified  as  a  species  from  Japan  (Hajek  et  al. 
1990).  Since  1989,  the  fungus  has  significantly  suppressed 
the  gypsy  moth  in  the  generally  infested  areas  of  North 
America.  Numerous  scenarios  have  been  proposed  for  how  it 
became  established  (Hajek  et  al.  1995),  but  we  probably  will 
never  know  for  sure. 

Fortunately,  E.  maimaiga  is  a  very  specific  natural 
enemy,  attacking  only  the  gypsy  moth  and  a  few  closely 
related  Lymantriidae.  Many  aspects  of  its  biology,  ecology, 
and  epidemiology  have  been  reviewed  by  Hajek  (1999).  The 
fiingus  passes  the  winter  as  spherical  resting  spores  in  forest 
soil.  In  spring,  at  about  the  time  when  gypsy  moth  larvae  are 
hatching,  and  if  the  soil  is  moist  enough,  some  spores 
germinate  to  produce  long  germ  tubes  that  form  pear-shaped 
primary  conidia  at  the  tip.  These  conidia  are  dispersed  in  the 
air  and  if  they  contact  a  gypsy  moth  larva,  will  break 
through  the  skin  to  infect  the  caterpillar. 

The  fungus  develops  rapidly  depending  on  temperature, 
eventually  killing  the  larva.  Then  if  relative  humidity  is 
close  to  100%,  the  fungus  will  sporulate  through  the 
caterpillar's  integument  and  produce  thousands  of  secondary 
conidia,  which  are  very  similar  to  primary  conidia.  These 
disperse  in  the  air  to  infect  other  larvae.  The  entire  process 
of  sporulation,  successful  air  transport,  and  infection 
depends  on  humidities  remaining  near  100%.  Near  the  end 
of  larval  development,  the  fungus  in  many  dead  infected 
larvae  form  thousands  of  resting  spores  rather  than  conidia. 
These  cadavers  then  weather  off  tree  trunks  and  release 
spores  into  the  top  soil  layers  where  they  over  winter. 
Weseloh  and  Andreadis  (2002)  have  shown  that  resting 
spores  may  remain  viable  in  soil  for  up  to  10  years,  thus 
providing  a  very  important  multi-year  reservoir  of  infection. 


Proper  weather  conditions  are  very  important  for  the 
fungus  to  be  effective.  Moist  soil  during  the  gypsy  moth 
larval  season  and  adequate  rainfall  at  the  times  of  conidial 
sporulation  must  occur  for  the  fungus  to  infect  larvae.  The 
number  of  resting  spores  in  the  soil  (resting  spore  load)  and 
the  number  of  gypsy  moths  available  for  infection  probably 
also  influence  fungus  activity.  Using  these  inputs,  colleagues 
and  I  have  developed  a  computer  model  that  simulates  the 
infection,  growth,  and  death  of  the  host  and  pathogen  over  a 
single  season  (Weseloh,  et  al.  1993,  Weseloh  1999,  2002). 
This  model  was  validated  with  independent  data  in  2002 
(Weseloh,  submitted).  Such  a  model  should  be  useful  in 
exploring  the  interaction  between  the  gypsy  moth  and 
E.  maimaiga,  and  might  even  be  helpftil  in  predicting 
seasonal  effectiveness  of  the  fiangus  if  some  indication  of 
future  weather  conditions  can  be  obtained.  To  encourage  its 
use  among  forest  managers,  I  have  developed  what  I  hope  is 
an  easily  used  version  of  the  model. 

This  publication  brings  together  information  about  the 
workings  of  the  model  and  the  kinds  of  inputs  necessary  for 
running  the  user-friendly  version.  The  latter  uses  less  exact 
estimates  for  input  data  than  does  the  validated  model,  so  I 
also  compare  outputs  of  the  two  versions.  The  model  is 
available  on  the  Internet  (at  http://www.caes.state.ct.us) 
and  can  be  downloaded  and  run  with  example  inputs  or 
your  own  inputs.  A  source  file  in  Borland  Turbo  C++  is 
also  provided  for  anyone  who  wishes  to  make  modifications 
to  the  program. 

Structure  of  the  Model.  The  model  was  developed  over  a 
number  of  years  and  reported  in  a  series  of  publications 
(Weseloh,  et  al.  1993,  Weseloh  1999,  2002).  These 
publications  should  be  consulted  for  a  description  of  model 
development  and  justification  of  equations  and  parameter 
values  used.  Only  the  current  model  will  be  described  here. 
The  model  simulates  the  interaction  between  E.  maimaiga 
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and  the  gypsy  moth  during  one  season  for  the  period  of  time 
that  gypsy  moth  larvae  are  available.  The  time  step  is  one 
day.  Larvae  are  assumed  to  hatch  over  a  period  of  16  days, 
starting  at  an  assigned  first  hatch  date,  such  that  the  number 
hatching  each  subsequent  day  varies  according  to  a  normal 
distribution  with  mean  of  7.5  days  after  the  start  date  and 
standard  deviation  of  2.5  days.  The  larvae  that  hatch  during 
a  given  day  are  assumed  to  be  one  cohort,  and  all  develop  at 
the  same  rate.  Cohorts  develop  according  to  a  degree-day 
model  with  threshold  temperature  of  12  C  (54  F).  They 
change  instars  at  the  following  degree-days  (dd)  after  hatch: 
(1-2  instar,  35  dd;  2-3  instar,  70  dd;  3-4  instar,  105  dd; 
4-5  instar,  150  dd;  >  5  instar,  240  dd). 

Resting  spores  germinate  when  the  soil  is  moist  enough. 
The  probability  of  infection  of  one  larva  in  a  cohort  by 
resting  spores  (Pr)  on  a  day  is: 

Pr  =  0.0035  R  [sin(  0.01431  M)]" 
where  R  is  the  percent  of  larvae  infected  in  a  soil  bioassay 
(see  inputs,  below),  and  M  is  a  forest  moisture  index  defined 
as  the  total  amount  of  rain  accumulated  over  the  present  day 
and  the  last  4  days,  with  no  more  than  10  mm  (0.4  inches)  of 
rain  allowed  per  day.  At  the  beginning  of  the  season  before 
secondary  conidia  are  produced,  the  number  of  larvae 
infected  from  a  given  cohort  on  a  given  date  is  determined 
by  multiplying  P,.  by  the  number  of  larvae  in  the  cohort. 
These  infected  larvae  then  form  a  new  infected  cohort  in 
which  the  fungus  develops  at  a  rate  (/•)  according  to 
temperature  (T)  between  7.5  to  26  C  (45-79  F): 

/•  = -0.162 +  0.032  7- 0.000752  T  " 

This  equafion  was  estimated  from  data  by  Hajek,  et  al. 
(1993).  From  the  rates  for  each  day,  the  model  determines 
the  date  when  larvae  in  the  infected  cohort  die. 

When  an  infected  cohort  dies,  the  number  of  conidia 
produced  per  m'^  of  forest  space  (i.e.,  the  volume  from  top  of 
canopy  to  ground)  per  day  (c)  varies  with  the  number  of 
larvae  in  the  cohort  (A';),  the  cohort  instar,  the  proportion  of 
the  day  for  which  humidity  is  above  98%  (//),  and  the 
number  of  days  after  death  (d): 

Instar  2:  c=  11639  Af,//exp(-0.36  cO 

Instar  3 :  c  =  24542  /V,  H  exp(-0.763  d) 

Instar  4:  c  =  32832  Af,//exp(-0.491  c/) 

Instar  5+:  c  =  204875  A',//exp(-0.859  if) 

If  c^  is  greater  than  8  days,  then  c  =  0.  Instars  5  and  6  are 
treated  as  one  instar  (5+)  because  the  equation  for  6th  instars 
was  similar  to  the  one  for  5th  instars  and  only  female 
caterpillars  generally  have  6  instars. 

Infection  of  a  cohort  by  secondary  conidia  depends  on 
conidial  concentration  (c)  as  determined  above,  and  a  size 
factor  for  the  larvae  in  a  cohort  (5) 

The  size  factor  is  determined  as  a  function  of  the  degree- 
days  of  development  (D)  of  the  healthy  cohort  fi-om  hatch: 
s  =  D-/(91613.9  +  0.6798D') 


Then  the  probability,  P,.,  that  one  larva  in  a  cohort  will 
become  infected  when  there  are  c  conidia  in  a  m^  volume  is: 

P,=  l-(  1  -0.0025)" 
where  s  is  the  above-mentioned  size  factor. 

The  total  probability  of  infection  (P/)  from  both  resting 
spores  (Pr)  and  conidia  (P^)  for  a  larva  in  a  cohort  is  then: 

Pr=Pr  +  Pc-P,Pc 
because  the  two  probabilities  are  independent.  The  model 
multiplies  this  total  probability  by  the  number  of  healthy 
caterpillars  in  a  cohort  to  give  the  number  of  larvae  in  a  new, 
infected  cohort. 

Model  Inputs.  The  model  inputs  are  daily  maximum  and 
minimum  temperatures  in  degrees  Fahrenheit  and  daily 
rainfall  in  inches  and  hundredths  of  inches  from  the  time  that 
gypsy  moth  larvae  begin  to  hatch  to  the  time  when  about 
half  of  them  have  pupated.  In  Connecticut,  this  is  from  about 
May  1  to  July  1 .  These  data  can  be  obtained  from  a  nearby 
weather  station,  but  because  rainfall  is  often  spatially 
localized,  it  is  best  if  daily  rainfall  is  measured  in  the 
location  of  interest.  These  data  are  adequate  for  determining 
degree-days  of  development  using  the  modified  sine-wave 
method  of  Allen  (1976)  for  interpolating  temperatures  within 
a  day  from  daily  maximum  and  minimum  values.  The 
proportion  of  a  particular  day  having  humidities  above  98% 
(//),  which  is  needed  to  determine  production  rate  of  conidia, 
is  estimated  from  daily  rainfall  on  the  present  day  (r,,  given 
in  mm),  and  the  last  3  days  {rj,  ri,  /■;)>  and  the  maximum  and 
minimum  temperatures  for  the  present  day  (in  degrees  C) 
(7;.,  T„,  respectively)  as://=  0.322998  +  0.1  sqrt  (  I. Or,,  +  0.5 
O  +  0.25  r,  +  0. 1 25  rj )  -  0.07  sqrt  (  7;,  -  r„  ) 

This  equation  was  obtained  after  a  good  deal  of  effort  to 
relate  rainfall  and  temperature  to  H,  and  the  fit  is  only 
approximate.  A  better  solution  is  to  monitor  actual  humidity 
levels  at  the  site  with  a  data  logger  capable  of  taking 
measurements  every  10  minutes  or  so.  The  validated  version 
of  the  program  uses  temperature,  rainfall,  and  humidity 
records  with  such  resolution  and  is  available  from  the  author 
(Ronald.  Weseloh@po. state. ct. us). 

Also  necessary  as  model  inputs  are  some  measure  of 
resting  spore  load  and  gypsy  moth  larval  density.  Resting 
spore  load  in  the  soil  is  best  measured  by  collecting  soil 
from  just  under  the  litter  layer  near  the  trunks  of  gypsy  moth 
favored  food  trees,  such  as  oaks,  mixing  these  together  for 
any  one  plot,  and  carrying  back  a  sample  to  the  laboratory. 
Forty  ml  sub-samples  of  the  soil  are  placed  in  a  2.5  x  14.0 
cm  plastic  Petri  dish  with  five  disks  of  12.5  cm  diam  filter 
paper  in  the  bottom  and  drenched  with  distilled  water.  After 
holding  for  1  week  at  17  C  and  15  hr  light  per  day,  25  2nd  or 
3rd  instar  gypsy  moths  are  placed  on  top  of  a  screen  disk 
placed  over  the  soil  and  held  with  the  soil  for  24  hr.  The 
larvae  are  then  reared  individually  in  40  ml  clear  plastic  cups 
filled  'A  full  with  gypsy  moth  diet  (Southland  Products,  Lake 
Village,  Arkansas)  and  held  for  2  weeks.  The  percentage  that 
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Table  1.  Comparison  of  seasonal  mortality  between  the  detailed,  validated  fungus  model  and  the  current,  user-  friendly 


version  for  various  areas  sampled  in  2001  and  2002. 


Seasonal  Mortality  prediction; 


Egg  masses/ 

Avg.  years  since 

Validated 

Current 

Year 

Location 

No.  Plots 

acre 

defoliation 

model  % 

model  % 

2001 

Harwinton,  CT 

2 

226 

>10 

99.2 

9 

2001 

Rhode  Island 

2 

2343 

1 

92.0 

100 

2001 

Mansfield  Hollow,  CT 

6 

2143 

5.9 

97.7 

98 

2002 

Mansfield  Hollow,  CT 

7 

261 

5.2 

62.3 

68 

die  from  fungus  is  a  measure  of  resting  spore  load  in  the  soil 
and  is  used  as  a  model  input. 

The  above  procedure  is  too  involved  for  common  use,  so 
the  current  model  version  uses  a  surrogate  measure  of 
resting  spore  load  based  on  the  number  of  years  since 
defoliation  last  occurred.  Weseloh  and  Andreadis  (2002) 
determined  that  resting  spore  viability  declines  linearly  over 
10  years.  Thus,  the  model  assumes  that  resting  spore  load 
can  be  expressed  as  the  greater  of  1 .0%  or  55.0%  -  5  y,  when 
y  is  the  number  of  years  since  the  last  defoliation  occurred  at 
the  site  of  interest. 

The  model  input  used  for  gypsy  moth  density  in  the 
validated  model  is  the  average  of  the  number  of  gypsy  moth 
larvae  found  on  25,  30-cm  (12  inches)  terminal  branches  of 
favored  food  plants  (such  as  oaks  or  witch  hazel)  that  can  be 
reached  from  the  ground.  This  count  is  made  about  3  weeks 
after  larvae  have  hatched.  The  current  model  assumes  that 
such  terminal  counts  can  be  derived  from  egg  mass 
counts/acre  by  multiplying  by  0.0007512  (Weseloh  2003). 
Thus,  the  input  to  the  model  is  the  number  of  egg 
masses/acre. 

In  the  course  of  developing  and  validating  the  model,  1 
found  that  local  dispersal  of  secondary  conidia  must  be  taken 
into  account  if  output  is  to  be  accurate.  Specifically,  the 
validated  model  works  much  better  if  the  number  of  conidia 
produced  each  day  from  several  plots  within  about  a  2  km 
(1 .2  miles)  radius  are  averaged  before  being  used  to  infect 
healthy  cohorts  (Weseloh  2002).  Within  such  a  small  area,  it 
is  likely  that  weather  conditions,  even  rainfall,  will  be 
homogeneous  enough  for  model  purposes.  But  resting  spore 
load  in  the  soil  and  gypsy  moth  larval  density  can  vary  by 
orders  of  magnitude  over  these  distances  (personal 
observations).  Thus,  it  is  necessary  to  obtain  defoliation 
history  and  egg  mass  density  for  five  to  ten  plots  within 
about  a  mile  of  the  location  of  interest.  The  values  are  then 
averaged  and  used  as  inputs  to  the  model. 

Model  Outputs.  The  model  generates  two  outputs:  (1)  the 
percent  of  living  larvae  infected  with  the  fungus  during  each 
day  of  the  season,  and  (2)  the  proportion  of  hatched  larvae 
that  are  healthy  on  each  day.  The  first  output  can  be  directly 
compared  to  actual  percent  infection  prevalence  in  the  forest 
as  sampled  by  collecting  and  rearing  larvae.  The  second 
output  is  most  useful  when  trying  to  evaluate  the  impact  of 


the  fungus  on  a  gypsy  moth  population.  The  program  always 
plots  these  values  in  a  graph  when  an  analysis  is  done  and 
displays  the  overall  seasonal  fungus-induced  gypsy  moth 
mortality.  If  you  name  an  output  file,  the  outputs  for  each 
day  are  saved  to  a  text  file  that  can  be  imported  into  a 
statistical  program  or  spreadsheet  for  further  analysis. 

Comparison  of  Model  Results.  In  2001  and  2002, 1  had 
available  high-quality  input  data  from  data-loggers  in  plots 
(Hobo  Pro  Rh/Temp,  Onset  Computer  Corp,  Pocasset,  MA) 
that  recorded  temperature  and  relative  humidity  every 
10  minutes,  and  a  data-logging  rain  gauge  (Onset)  that 
recorded  precipitation.  Also,  resting  spore  loads  from 
bioassay  results  and  terminal  counts  were  available  both 
years.  Thus,  I  had  data  suitable  for  using  the  validated  model 
and  from  it  I  obtained  values  of  average  seasonal  mortality 
in  four  different  areas,  three  in  2001  and  one  in  2002.  These 
estimates  are  compared  to  estimates  for  the  same  areas  using 
the  current,  user-friendly  version  (Table  1). 

Only  for  the  Harwinton  area  in  200 1  were  the  results  for 
the  two  models  substantially  different.  Part  of  the  reason  is 
that  actual  resting  spore  load  was  much  higher  in  this  area 
than  was  estimated  by  defoliation  history.  However,  even 
when  I  assumed  defoliation  had  occurred  the  year  before, 
seasonal  mortality  in  Harwinton  for  the  current  model  only 
increased  to  58%.  The  Harwinton  result  shows  the 
inadvisability  of  using  data  from  only  a  few  locations  (in  this 
case,  two). 

Future  Uses.  The  current  version  of  the  model  should  be 
relatively  easy  to  use  because  the  required  biological  input 
data  (defoliation  and  egg  mass  records)  are  often  readily 
available.  Fortunately,  the  user-friendly  version  often 
approaches  the  accuracy  of  the  validated  version.  However, 
weather  information  is  also  needed  for  input,  and  because 
weather  cannot  be  predicted  in  any  precise  way,  neither  can 
future  fiingus  effects.  But  the  model  can  be  used  to  explore 
potential  fungal-host  population  dynamics  under  a  variety  of 
weather  conditions.  This  could  be  of  some  value,  because  if 
resting  spore  load  and  gypsy  moth  densities  are  high,  the 
fungus  is  most  likely  to  control  the  gypsy  moth  unless  the 
season  is  expected  to  be  extremely  dry  (Weseloh  2002). 
Various  weather  scenarios  can  be  explored  to  gain 
information  about  what  is  likely  to  happen.  The  model 
should  also  be  useful  in  the  further  exploration  of  gypsy 
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moth-flingus  population  ecology.  The  prevalence  of 
infection  output  can  be  used  to  check  that  model  output  is 
similar  to  actual  fungal  prevalence,  while  the  survival  output 
can  be  used  to  evaluate  fungus  effectiveness.  Thus,  survival 
data  can  be  derived  from  disease  prevalence  data,  something 
that  is  not  usually  possible.  This  ability  makes  the  model 
output  valuable  for  studying  the  population  dynamics  of  this 
system. 

Availability.  The  program  runs  on  Microsoft  Windows© 
and  can  be  downloaded  from  the  Connecticut  Agricultural 
Experiment  Station  Website  (http://www.caes.state.ct.us). 
The  downloaded  file  (maimaiga.zip)  can  be  uncompressed 
using  WINzip  or  PKzip.  For  best  results,  place  the  expanded 
files  in  a  separate  directory  called  \fungus\  so  the  example 
weather  files  can  be  accessed  easily.  Among  the  files  is 
"readme.txt",  that  has  detailed  information  about  the  model. 
To  run  the  model,  start  the  file  "fungus  1  .exe".  Alternately,  a 
disk  with  all  relevent  files  and  directions  for  installations  is 
available  from  the  author. 
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